AD  6772  72 


TG— 1007 
AUGUST  1968 
Copy  No.  2 


Technical  Memorandum 

APPROXIMATE  FUNCTIOHS  AS 
PARTICULAR  SOLUTIONS  IN 
THERMAL-STRESS  ANALYSIS  OF 
AN  OGIVAL  RADOME 

by  MANFORD  B.  TATE 


i 

i 


THE  JOHNS  HOPKINS  UNIVERSITY  •  APPLIED  PHYSICS  LABORATORY 


Thu  document  hoi  b»n  opproxtd  lot  public 
r*l«at«  and  tale,  iti  diitribution  it  wnlimittd. 


TG— 1007 

AUGUST  1968 


Technical  Memorandum 

APPROXIMATE  FUNCTIONS  AS 
PARTICULAR  SOLUTIONS  IN 
THERMAL-STRESS  ANALYSIS  OF 
AN  OGIVAL  RADOME 

by  MANFORD  B.  TATE 


THE  JOHNS  HOPKINS  UNIVERSITY  •  APPLIED  PHYSICS  LABORATORY 
8621  Georgia  Avenue,  Silver  Spring,  Maryland  20910 

Operating  under  Contract  NOw  62*0604-c,  Bureau  of  Naval  Weapons,  Department  of  the  Navy 


This  document  hot  bssn  opprovsd  for  public 
rcWatff  and  lolt;  its  distribution  it  unlimited. 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORT 
Silvch  Shmno  marylano 


ABSTRACT 


Functions  of  approximation  to  particular  solutions  that  occur  in  the 
ogival  radome  thermal-stress  problem  are  presented.  It  was  found  by  numerical 
comparisons  with  the  generating  differential  equations  that  the  approximate 
closed-form  solutions  display  an  error  range  of  no  more  than  plus  or  minus 
one-half  percent. 

The  solutions  were  derived  to  gain  two  major  advantages:  first,  a  reduction 
of  analytical  complexity  lessens  the  chance  of  computational  error;  and,  second, 
certain  unwieldiness  in  numerical  work  is  reduced  or  eliminated. 

Computer  results  are  tabulated  for  repeated  future  application  in  the 
evaluation  of  thermal  stresses  in  blunt  and  pointed  radomes  of  compound-ogive 
configuration. 
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I .  INTRODUCTION 


General  analyses  of  shell  stresses  are  developed  in  textbooks,  References 

(I)  to  (5),  inclusive,  and  Boley  and  Weiner  (4)  discussed  thermal  stresses. 

For  heat-variant  material  properties  in  the  wall  of  an  ogival  radome,  axisym- 
metric  thermal  stresses  were  studied  in  Reference  (6). 

In  References  (7)  and  (8),  Rivello  took  up  the  problem  of  thermal  stress 
in  cylindrical  sandwich  shells,  and  Dailey  (9)  employed  the  stiffness-matrix 
method  for  axisymmetric  shells  of  revolution.  Weckesser,  Hallendorff  and 
Suess  (10)  investigated  materials  for  radome  construction  in  high-temperature 
applications.  Thick-walled  conical  shells  were  examined  by  Weiss  in  Reference 

(II) . 


Pyroceram  9606  test  data,  temperatures  around  the  nose  of  a  blunt  radome, 
the  Von  Karman  shape,  temperature  distributions  in  a  test  radome,  and  compound- 
ogive  profiles  were  investigated  in  References  (12)  to  (16),  inclusive,  in 
association  with  the  present  subject. 

From  numerical  studies,  it  was  found  that  particular  solutions  based  on 
applied  temperature  distributions  obtained  as  infinite  series  made  computer¬ 
ization  more  complex,  inhibited  analytical  continuity,  and  increased  an  amount 
oi  computation  that  was  already  rather  cumbersome.  To  offset  these  difficulties, 
closed-form  approximate  solutions  applicable  to  particular  cases  are  developed 
in  the  present  text.  Moreover,  computerized  data  on  these  solutions  are 
tabulated  for  future  reference. 


The  ogival  radome  wall  profile  is  sketched  on  Figure  1,  where  dimensions 
and  coordinate  variables  are  defined  also.  A  differential  element  of  the 
wall  is  drawn  in  Figure  2.  On  it  are  shown  normal-stress  resultants 
bending  moments  (M^ ,  M^),  and  the  shearing-stress  resultant  (Q^). 
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Fig.  2  RADOME-WALL  ELEMENT  WITH  STRESS  RESULTANTS  AND  BENDING  MOMENTS 
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III.  NOMENCLATURE 


A,B:  Constants 

E:  Young's  modulus  of  elasticity,  psi 

J,K:  Constants 

L:  Operator 

M:  Bending -moment  resultant,  ippi 

N:  Normal-stress  resultant,  ppi 

Q:  Shearing-stress  resultant,  ppi 

R:  Radius,  inches 

V:  Wall-slope  function,  radians  per  inch 

Z:  Geometric  axis  of  radome 

c:  Radome-wall  half  thickness  (c  ■  h/2),  inches 

h:  Radome-wall  thickness,  inches 

t:  Span  length  of  radome,  inches 

q;  Auxiliary  shear  function,  ppi 

r:  Radius,  inches 

A:  Denominator 

0:  Wall-bending  parameter,  radians 

8:  Coordinate  angle  of  rotation,  degrees  or  radians 

v:  Poisson's  ratio 

X:  Radome-wall  curvature  change,  radians  per  inch 

Coordinate  angle  of  azimuyh,  degrees  or  radians 

The  following  notations  are  employed  as  subscripts: 

a:  Anterior  or  outer  surface  of  radome 

b:  Base  of  radome 

c:  Central  surface  of  radome  wall 

g:  General 

i,j:  Numerical  indices:  0,  1,  2,  3,... 

o:  Origin  or  initial  value  (zero) 

p:  Particular 

s:  Secondary  or  inner  surface  of  radome 

v:  Pertaining  to  wall  slope  (V) 

R,6,t:  Spherical-coordinate  directions 


4 


THE  JOHNS  MOW  INS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
Suvn  iwino  ma*vi_amo 


IV.  REFERENCE  EQUATIONS 


The  differential  equations  derived  in  Reference  (6)  are: 

r  K  R 

(La+Y)(q-H^)  =  EhV  +  (l-v)^  +  8 

c  rc  ” 


(La -v) (EhV) 


(v2-^4) 


+ 


R 

c 


q 


(1) 

(2) 


In  these  expressions,  q  is  the  auxiliary  shear  function  defined  by  equation  (3) 
wherein  (fy  is  resultant  or  total  shear,  V  represents  change  in  slope  of  the 
radome  wall,  v  is  Poisson's  ratio,  E  is  Young's  modulus  of  elasticity,  L  is  the 
operator  defined  below  in  equation  (4),  rc  and  Rc  are  the  radii  pictured  in 
Tigure  1  where  the  coordinate  angle  *  also  is  shown,  and  the  other  terms  are 
constants . 


!s3t 

Rc  sin* 


(3) 


L‘q  = 


R^  sin*  d* 


j2  .  qR  cosi|i 

&  +  cot,  £  C 


d*  r£  tart* 


(4) 


Solutions  of  differential  equations  (1)  and  (2)  are  obtained  as 


q 


q  + 
4p 


V 


V  =  V  +  V 
p  g 


(5) 


where  subscripts  p,  g  denote  the  particular  and  general  solutions,  respectively. 
Our  purpose  herein  is  to  develop  functions  of  approximation  to  the  particular 
solutions . 


V.  TOTAL  SHEAR 


The  total  shear  (Qi )  can  be  calculated  with  equation  (3)  from  function  q 

_L  J  _  _  JT 1  _  1^  J 1  _  _  1 _ _ _ 4 1 1 1 


for  which  the  particular  solution  qp  is  given  approximately  by 


%  =  V 


rc^B  Rc0*y  +  ^cos*) 
R  +  r  tan* 


(6) 


wherein  the  are  constants,  and  the  remaining  quantities  are  illustrated  in 
Figures  1  and  2. 
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VI.  WALL-SLOPE  CHANGE 


the  particular  part  V  of  slope  change  V  of  the  radome-wall  profile  is 
calculated  approximately  with 


V 

P 


rcV  Rc(^,+  **  cost) 

R  r  tank 

c  c  ^ 


(7) 


where  the  B!  are  constants,  and  the  radii  (r  ,  Rc)  and  the  coordinate  angle  are 
identified  In  Figure  1. 


VII.  NORMAL-STRESS  RESULTANTS 


The  normal -stress  resultants  are  expressed  in  terms  of  the  foregoing  relations 
as  follows: 


q  R  cost 


(8) 


N 


0P 


B  R  R  r 

_  ■  c  .  c 

h  -  -7-  cost  -  —  C 

.  c  J  c 


CSCSt  + 


Rc  cost 

rc  tan* 


(1^+  ^  cost) 


(9) 


which  contain  the  constants  and  variables  referred  to  in  connection  with  equation 

(6). 


VIII.  BENDING  CURVATURES 


The  components  of  curvature  change  produced  by  bending  of  the  radome  wall 

are : 


V  cost 

y  «  _E _ 1 

*0p  r 

c 


(10) 


\  v‘ 

cost  + 

R  cost 

A  P  T 

r  "  R 

csc2t  +  - — —— 

r  tartf 

c  c 

—  - 

c 

c 

(*;+  y  cost) 


(11) 


which  contain  the  constants  and  variables  referred  to  in  connection  with  equation 
(7). 
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IX.  TEMPERATUP  J  FUNCTIONS 


As  developed  in  Reference  (6),  temperature  functions  are  defined  by  the 
following  expressions. 


at  =  Eet/(l-v) 


(12) 


(13) 


(14) 


With  x  =  y/c,  the  thermal  strain  («t)  obtained  in  Reference  (12)  can  be  written 
as  follows. 


6t  =  (-241.2  +  3.445T)  x  10"6 ,  *1  S  x  S  &  (15) 

=  (-119.1  +  2.958T)  x  1(T6,  0  £  x  £  0.5  (16) 

=  (+190.3  +  2 . 22 IT)  x  10"5 ,  0.5  £  x  £  1  (17) 


The  preceding  relations  are  applicable  to  the  entire  radome  for  the  temperature 
range  of 


70°F  £  T  £  1400°F 


(18) 


and  are  used  together  with  T,  which  is  expressed  by 


T  =  T(x,*)  =  f(x)  T  (*) 

a 


(19) 
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where  is  the  outer-surface  temperature  distribution  plotted  on  Figure  3.  The 
spanwise  linear  temperature  distribution  shown  in  Figure  3  as 


Ta  "  Ta^>  "  589  +  705  cos*  (20) 

was  employed  in  the  reported  computations.  And  f(x)  giving  the  distributive 
character  over  the  thickness  of  the  wall  is 


3 

f (x)  =  V  fnxn  *  0.257  +  0. 358x  +  0.305X2  +  0.084X3 

n*b 


(21) 


wherein  the  polynomial  coefficients  f  are  shown  in  the  expanded  form  on  the 
right. 

When  the  integrations  of  equations  (13)  and  (14)  are  carried  out  with 
equations  (12)  and  (15)  to  (21),  the  resultant  temperature  functions  are  obtained 
as  written  below. 


R 


Nte  “  Ko  +  * 

cost  , 

Ntt  -  Nte  +  f 

T  c 

(Kg+  %  cost  ) 

(22) 

R 

Mte  *  Jo  +  Jx 

cost  , 

Mtt  "  Mte  +  T 

r 

(Ja+  Jg  cost) 

(23) 

c 


The  constants  J.,  are  the  same  as  those  appearing  in  equations  (1)  and  (2) 


and,  by  the  methods  just  described,  were  'found  to  be 

Kq  =  4.648Eh/(l-v)  x  id*  ,  Kj  -  6 . 792Eh/(l-v)  x  10*  (24) 

Kg  «  4.393Eh/(l-v)  x  107  ,  K,  =  3.420Eh/(l-v)  x  l(f  (25) 

JQ  =  1.269Eha/(l-v)  x  lCf,  -  0.989Eha  / (1-v)  x  ltf  (26) 

Ja  =  1.704 Eh8  / (l-v)x  107  ,  ^  -  2 . 195Eha  /(1-v)  x  107  (27) 
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and  the  material  properties  were  evaluated  from  test  data  on  Pyroceram  9606  in 
Reference  (12);  i.e., 


E  *  16,400,000  psi  v  *  0.244 


(28) 


where  E  is  accurate  within  four  percent  and  v  within  2.5  percent  over  temperature 
range  (18). 


X.  DETERMINATION  OF  CONSTANTS 


By  putting  expressions  (6)  and  (7)  into  (1)  and  (2),  we  secure  formulas 
for  numerical  determination  of  constants  and  B^.  These  constants  are  com¬ 
puted  in  the  following  sequence. 


J,  (l-v2)(Jl+K1Rc) 
^  "  Rc  (14-6*  )  Rc 


(l*v)  (v3 40*  )  (Jj  -HCj  R£) 
(140*  )Eh  Rc 


(29) 


(v24B4)(J34KjR  ) 

*  — S3 - • 


\  “  Ag-K,  (30) 


0*  -2 


4 


2  £ 

B4  -2  4  ’ 


0*  -  -Vs  4  12(l-va)(Rc/h)a  (31) 


r  (3  4  2  cot3*) 
(3  4  cot3*)  4  -r-.IS* - 


+  4s 

^-2  4  t» 


V 


-'mean 

^  -  (A,+  v)^ 

Eh 


vAb  4  2% 

EL*  *  —  ' 

^  Eh 


(Vv)Ka 

4+  e4 


(v3 4  04)Kg 
^  "  (4|  4  0*  )Eh 


1  4 


R  sirt3* 
c 


mean 


(32) 

(33) 

(34) 

(35) 

(36) 
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For  slender  ogives,  cotf  and  rc/Rc  are  small  and  sin|  approaches  unity. 

The  terms  (32)  and  (36),  therefore,  have  slight  variation,  and  one  can  employ 
mean  values  computed  as  arithmetical  averages  with  initial  and  final  values 
of  the  coordinates.  For  the  bicentric-ogive  radome  described  in  Rfference  (16), 
we  have 


rcl 

0.238", 

=  64°  19' 

23" 

(37) 

*cb 

6.625", 

♦b  =  9°° 

(38) 

= 

0.25", 

R  =  64.679" 
c 

(39) 

and  the  numerical 

results  are 

listed  in  Tables  1  and 

2. 

XI.  COMPUTER  RESULTS 


As  part  of  the  overall  computer  program,  values  of  the  functions  occurring 
in  the  particular  solutions  were  calculated  in  Program  Problem  164.  The 
associated  constants  are  reported  in  Table  1;  and  functional  numerical  data, 
in  Table  2.  All  of  the  computer  results  that  are  reported  in  Tables  1  and  2 
were  computed  and  stored  in  double  precision  but  shown  here  in  single  precision. 
The  double  precision  numbers  are  needed  in  simultaneous  solutions  for  general 
constants  of  integration,  which  are  evaluated  from  specified  boundary  conditions 
imposed  at  the  juncture  between  nosecap  and  mainbody  and  at  the  base  of  the 
radome  where  it  is  joined  to  a  missile  body. 

XII.  DISCUSSION 

Formulation  of  closed-form  approximate  particular  functions  (qp,  Vp,  Ng  , 
Xi_)  are  presented  as  equations  (6),  (7),  (9),  and  (11)  for  the  mainbody  of 
tnl  bicentric-ogive  radome  described  in  Reference  (16).  They  were  programmed 
for  numerical  evaluation  at  thirteen  points  along  the  length  of  the  mainbody 
lying  in  the  interval  $  £  90°,  where  *  64°  19'  23"  and  the  base  of 

the  radome  is  located  at  90°  (Figure  1). 

The  results  are  reported  for  future  reference,  because  they  are  used  in 
each  analysis  of  thermal  stress  and  deformation  irrespective  of  the  boundary 
conditions  imposed  on  the  radome.  Specific  values  occur  in  the  fulfillment 
of  boundary  requirements  which  provide  simultaneous  equations  that  occur  in 
calculation  of  integration  constants  appearing  in  the  general  integrals  of 
differential  equations  (1)  and  (2). 

One  can  observe  from  Table  2  that  qp  increases  monotonical ly  from  the 
nosecap  junction  (at  ^ )  to  the  radome  base  (at  90°).  But  total  shear  in 
the  particular  solution  behaves  oppositely  as  computed  with  the  following 
relation. 
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Table  1.  Constants  in  Closed-Form  Particular  Solutions 


Symbol 

Value 

Symbol  Value 

e4 

755,387.78 

% 

0.460,533  x  10-s 

A* 

2.24 

Bb 

2.067,534 

4* 

3.25 

h 

0.629,548  x  10“6 

755,385.79 

Bb 

0.799,994  x  10“6 

B4 

+  taa 

755,392.81 

Bb' 

0.110,662  x  10-6 

B* 

-2  +  4 

755,396.35 

* 

6.795,812 

a; 

0.581,105  x  10-6 

*B 

1.859,424 

V 

0.453,511  x  10-6 

Table  2.  Functional 

Values  of  Particular  Solutions 

j 

t 

qp  x  1(f  * 

N0P  X  1Cf 

V  x  10*  R  x*  x  id4 

(d-m-s) 

Eh 

Eh 

1 

64-19-23 

0.328,962,85 

18.083,065  -1 

.040,084,8  119,874,77 

2 

65-06-39 

0.605,555,02 

20.680.353  -0 

.438,696,15  14,835,188 

3 

66-45-43 

1.191,282,1 

19.793,691 

.296,856,92  0.662,232,44 

4 

68-31-56 

1.781,837,4 

18.415,570 

.312,404,94  -  1.191,796,8 

5 

70-27-15 

2.373,419,8 

16,851,096 

.361,440,69  -  1.621,952,3 

6 

72-01-27 

2.817,373,3 

15,549,685 

.406,894,98  -  1.671,379,0 

7 

73-44-23 

3.2^1,446,1 

14.110,889 

.456,230,82  -  1.611,040,9 

8 

75-27-58 

3.664,639,9 

12.648,412 

.503,095,07  -  1.492,807,3 

9 

77-25-18 

4.067,873,5 

10.976,944 

.551,223,87  -  1.322,628,6 

10 

79-06-45 

4.370,315,4 

9.521,078,6  - 

.587,846,96  -  1.157,576,3 

11 

81-06-53 

4.672,765,4 

7.785,791,5  - 

.624,688,70  -  0.948,900,12 

12 

83-43-13 

4.975,213,2 

5.513,499,7  - 

.661,404,10  -  .663,819,83 

13 

90° 

5.277,579,8 

0.001/19,1  - 

.694,981,63  +  .056,732,43 
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qpRc  slt*» 

r 

c 


Here,  the  sine  increases  from  sin^fj  =  0.901250  to  sin  *.  ■  1;  and  r£ 
0.003679  to  0.102429.  Therefore 


=  80.586 ,510Eh  x  10-*  =  3.304  ppi 

^ p(90°)  =  51 .524 ,273Eh  x  10-*  =  2.112  ppi 
and  similar  comparisons  can  be  made  for  the  other  functions  in  Table  2 
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XIII.  CONCLUDING  REMARKS 


Approximate  functions  for  particular  solutions  of  governing  differential 
equations  are  presented  in  the  text.  They  were  developed  to  avoid  cumbersome 
infinite  series  and  to  express  the  relevant  functions  in  closed  form.  The 
generating  equations  are  satisfied  (numerically)  by  the  derived  expressions 
such  that  discrepancies  lie  in  a  range  of  plus  or  minus  one-half  percent. 

Major  advantages  of  these  solutions  are  simplification  that  reduces 
the  likelihood  of  numerical  errors  and  certain  unwieldy  characteristics  in 
numerical  work  are  minimized  or  eliminated. 

Computerized  data  on  ogival  radome  thermal-stress  functions  ate  collected 
that  can  be  used  in  future  stress  calculations. 
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